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Abstract 

We study the homogenization of a Schrodinger equation in a pe- 
riodic medium with a time dependent potential. This is a model for 
semiconductors excited by an external electromagnetic wave. We prove 
that, for a suitable choice of oscillating (both in time and space) po- 
tential, one can partially transfer electrons from one Bloch band to 
another. This justifies the famous "Fermi golden rule" for the transi- 
tion probability between two such states which is at the basis of various 
optical properties of semiconductors. Our method is based on a com- 
bination of classical homogenization techniques (two-scale convergence 
and suitable oscillating test functions) and of Bloch waves theory. 

Dedicated to the memory of Frederic Poupaud. 



1 Introduction 

This work is devoted to the mathematical justification of a problem of mean 
field approximation in solid state physics. More precisely, we study the ho- 
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mogenization of the following Schrodinger equation 

S^-An,+ (e-2c(^)+4(t,x))«, = inM^x(0,T) 



Ue{t = 0,x) = u%x) 



in 



t>N 



where < T < +00 is a final time. The potential c{y) is a real, bounded 
and periodic function defined for y G (the unit torus). Equation ^ is 
the so-called one-electron model for describing the electrons in a crystal or 
in a semiconductor, characterized by the periodic potential c{y) [2|, (H], [25] . 
An exterior field is also applied to the sample: in first approximation it is 
described by another real potential de{t,x) which depends also on the time 
variable in contrast to c (see Section El for the case of an electromagnetic 
potential). In JT]) the size of the background microscopic potential is of order 
while that of the exterior macroscopic potential is of order e°, so is just 
a small perturbation (the relative size of which is the square of the period). 
It can also be seen as a control acting on the semiconductor and its specific 
dependence on e will vary with the initial condition and the desired final or 
target state at time T (see ^ below). For example, the exterior potential 
can be light illuminating the semiconductor: if its energy is high enough, it 
can excite electrons from the valence band to conduction band. This effect 
is called optical absorption. Its converse effect (light emitted by electrons 
going from the conduction band to the valence band) is at the root of many 
important devices such as lasers, light emitting diodes and photo-detectors 

[HI, m. m 

Remark that Planck's constant has been normalized to unity in (jlj. We 
emphasize that the e-scaling in JH) is not the usual semi-classical scaling 
for Schrodinger equation [H], [H], [12], [iHl, [H], [201 which would involve 
a coefficient in front of the time derivative. Instead, it is the scaling 
of homogenization as in [2], [2], [1]. In physical terms it corresponds to an 
asymptotic study for much longer times than in the semi-classical limit (we 
refer to [3] for a more complete discussion of the scaling). 

Let us describe a typical example of our results. We introduce the so- 
called Bloch or shifted cell problem. 



-(div^ + 2ine) ({Wy + 2^7r^)^„) + c(y)^„ = XnWn in T 



where 9 G is a parameter and {\n{9),tpn{y, 9)) is the n-th eigencouple. In 
physical terms, the Bloch frequency 9 is the quasi momentum and the range 
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of A„(6'), as 9 varies, is a Bloch or energy band. We consider an initial data 
which is a wave packet of the type 

w°(x)=V^„(^,^'^)e2-^t;°(x), (2) 

and we would like to attain a final state at a different momentum 9"^ and 
energy \m{,9'^) 



ui{x)=^Pm[-,9"^)e''''—v'ix). (3) 
For achieving this goal we choose an oscillating potential 

,- (A,n(6)'")-A„(9"))t ni^ (e"'-e")-^ \ , / . X 



d,{t, x) = ^ [e' e''^ ^ j d (t, X, ^ j , (4) 

where 3? denotes the real part, and d{t,x,y) is a real potential defined on 
[0, T] X X T'^. Formula (jH) for the potential is quite natural: the oscil- 
lating phase in time corresponds to the energy difference between the initial 
and final state while the oscillating phase in space corresponds to 

the momentum difference. In other words, the potential puts energy and mo- 
mentum in the system so as to have global conservation of these quantities. 

In truth we cannot reach (even approximately) the desired final state (jSJ . 
Instead we end up with a mixed state, combination of ^ and Q. Under 
the assumption that 6'" and 6*™ are critical points of their non-degenerate 
energy levels (which implies that the group velocities of the wave packets ^ 
and ^ vanish) and are non-resonant (which means that no other state are 
excited by the external potential), see ^ and ^ below, we shall prove in 
Theorem 13.21 that the solution of (ID satisfies 



u,{t,x) ^ e^^e2-^V^„ v4t,x)+^^^e''^^^Pm (^,^") vUt,x), 

where the macroscopic profile (f„,fm) is the unique solution of the following 
Schrodinger homogenized coupled system 

Ov 

- div {AlVvn) + <„(t, x)vra = ^ in X (0, T) 



Ot 

Ov. 



l- 



" - div {AlyUm) + Cn(^, x)Vn = () iu X (0, T) (5) 



Ot 



Vnit = 0,x) = f°(x) in 
^ Vmit = 0,x) = in 
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with homogenized coefficients A^,A^ and c?*^ = d^^. The tensors A* and 
are the inverses of the effective masses of the particles corresponding 
to the initial and desired final state. The coupling coefficient d*^ is given 
by a formula (see (f2n|l and Remark l3.3p known as "Fermi golden rule" 
P], [23]. The homogenized system (0) is a model for light absorption in 
semiconductors. When 6'" = 6'™ one talks about "direct" absorption, and 
when 6*" 7^ 6"" about "indirect" absorption [Z|, [HI, [23, [23- In truth, one 
does not find ((HI), as it stands, in the physical literature where instead a 
simpler semi-classical picture is used. Specifically, physicists talk about the 
transition probability between the two states ^ and 0, which is precisely 
equal to the squared modulus of the coupling coefficient d*^. 

To obtain the homogenized limit ^ we follow the method introduced 
in [2], [SI- The main idea is to use Bloch wave theory to build adequate 
oscillating test functions and pass to the limit using two-scale convergence 

H, cni 

The content of this paper is as follows. Section [21 is devoted to recalling 
basic facts about Bloch waves and two-scale convergence as well as stating 
our main assumption. Section El gives our main result of homogenization. 
Section HI focuses on a generalization of the Schrodinger equation ^ which 
takes into account an electromagnetic field. Finally Section |S1 is concerned 
with a resonant case where more than two states are coupled. 

2 Bloch spectrum and two-scale convergence 

In this section we recall some results on Bloch waves and two-scale con- 
vergence, and we introduce our main assumptions on the initial and target 
states. 

We assume that the potential c{y) is a real measurable bounded periodic 
function, i.e. belongs to L°°(T^), while the potential d{t,x,y), appearing 
in is real, measurable, uniformly bounded, periodic in y and smooth in 
(t, x). We recall that, for given 9, the Bloch (or shifted) spectral cell equation 

- (div^ + 2t7ce) ((Vy + 2z7r^)7/'„) + c{y)4jn = XnWn in T^, (6) 

admits a countable sequence of real increasing eigenvalues (A„)„>i (repeated 
with their multiplicity) and normalized eigenfunctions {ipn)n>i, with ||V'ra||L2(T'V) = 

I, since its Green operator is a compact self-adjoint complex- valued operator 
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on L^(T^). The dual parameter 9 is called the Bloch frequency or quasi 
momentum and it runs in the dual cell of T^, which, by our choice of nor- 
malization factor 2ti in the phase factor, is again the unit torus T^. In other 
words, by periodicity it is enough to consider 6 G T^. For more details on 
Bloch waves, see e.g. [101, [Hli [21 • 

In the sequel, we shall consider two energy levels n, m > 1 and Bloch 
parameters 6"-, 6"^ G such that the eigenvalues A„(6'") and Xm{6"^) satisfy 
the following assumption 



Of course, we assume that, either n m, oi n = m and 9"- ^ 9'"^ (if n = m 
and 9" = 9"^, then = d(t,x,x/e) and this case was already treated in 
[3J). The simplicity assumption, i.e. part (i) of is generic and simplifies 
considerably the analysis. In particular, it implies that Xp{9) is infinitely 
differentiable in a vicinity of 9^, and one can introduce the group velocity 
V6iAp(6'). The criticality assumption, i.e. part (m) of (jH), is physically relevant 
when states at the bottom or top of Bloch bands are considered. For a 
discussion of this type of assumptions, as well as possible weaker ones, we 
refer to |[H]. 

We also add a non-resonant assumption 



(m) for any p>l, \{29'^ - 9"") ^ 2A„(r) - A™(^™). (8) 



The interpretation of assumption ^ is the following. The oscillating poten- 
tial dg, defined by has been designed to transfer the initial state with 
(quasi) momentum 9"' and energy A„(6'") to the target state 9'^,\m{9'^). 
The only requirement is that momentum and energy are conserved during 
this process. Actually there is another possible state that can be reached un- 
der the conservative action of d^, namely the state with momentum 29"- — 9"^ 
and energy 2A„(6'") — \m{9"^). In order to simplify the analysis, assumption 
(jHI) forbids this additional state as a standing wave solution of {TJ without 
exterior potential. Section El explores the resonant case where ^ is not 



Under assumption it is well-known [T^ that one can make a choice of 
an eigenvector such that the n-th eigencouple of © is smooth in a neighbor- 
hood of 9^. Introducing the operator An{9) defined on L^(T^) by 




(i) \{9^) is a simple eigenvalue, 
{ii) 9^ is a critical point of \{9) i.e., Ve\{9^) = 0. 



(7) 



satisfied. 
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we differentiate ^ with respect to 6. Denoting by {ek)i<k<N tfie canonical 
basis of and by {9k)i<k<N the components of 9, the first derivative satisfies 

K{0)^ = 2t7rekiVy+2t7re)^n+{diYy+2tTTe) (2z7refc^„) + i^(^)V'n, (10) 
and the second derivative is 

A„(^)^^ = 2z7refc ■ (V, + 2m^)^ + (div, + 2ni9) (szTre^^ 



dOk ' ' d9i dOi ^ ' ddk 
(9^ A 



(11) 

Under assumption (0) we have V6iA„(6'") = 0, thus equations (|TTHl and (|TI| 
simplify for 9 = 9^ and we find 

where satisfies 

K{9nCn = ek ■ (V, + 2t7r9^)tlj^ + (div, + 2i7r9^) (e^^^) in T^, (13) 

and Xn satisfies 

Ki&nXn = efc ■ (V, + 2tn9^)C + (div, + 2^71^) (e,Ci) 

+e, ■ (V, + 2^n9-)C + (div, + 2^n9^) (e,C) ^^^^ 

We know that ^^id Xn exist since they are defined by (fT2|l as multiple 
of the derivatives of ipn with respect to 9 (and assumption (0) ensures that 
ipn is indeed differentiable) . However, if we forget for a moment definition 
(fT2| . the existence and uniqueness of the solutions to (fT3ll and (fT4ll is not 
at all obvious. Since the operator A„(6'") has a non empty kernel spanned 
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by ipn, one should apply the Fredholm alternative: equations (|T!1|l and (fn|l 
admit a unique solution (up to the addition of a multiple of ipn) if and 
only if their right hand side are orthogonal to ipn (i-e. satisfy the Fredholm 
compatibility condition). This compatibility condition is not immediately 
satisfied. Actually, it gives new informations which are a consequence of the 
previously established existence of and Xn- particular, the compatibility 
condition of yields a formula for the Hessian matrix VeV6iA„(6'") in terms 
of ipn and that we shall use later (see (f3T|) V 

Remark 2.1 All our results can be generalized if we replace the Laplacian 
in (Cp by the more general operator div(A(?/)V-) where A{y) is a symmetric, 
bounded, periodic and uniformly coercive matrix. In this case, the Bloch 
spectral cell problem ^ becomes 



A tensor A{y) ^ Id may be interpreted as a periodic metric. It makes sense 
for the study of wave propagation in a periodic media (see e.g. 12]). 

Finally we recall the notion of two-scale convergence introduced in p], 



Proposition 2.2 Let be a sequence uniformly bounded in L^(]R^). There 
exists a subsequence, still denoted by u^, and a limit UQ{x,y) G L^(M^ x T^) 
such that two-scale converges (weakly) to Uq in the sense that 



for all functions (p{x, y) G (M^; C#(T^)) . 

Notation: for any function (p{x,y) defined on x T^, we denote by 0^ 
the function 0(x, |). 

3 Main result 

Due to our assumptions on the coefficients, if the initial data u° belongs to 
H\R^), there exists a unique solution of the Schrodinger equation (QJ in 
C ([0,T]; H^(R^)) which satisfies the following a priori estimate. 



(divj, + 2i7re) {A{y)(yy + 2i7i9)ijJ + c{y)i:n = A„(^)^„ in T 



<N 
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\u^{t, x)\'^dx = 0. 



Lemma 3.1 There exists a constant C > 0, which depends on T but not on 
e, such that the solution of (CJ) satisfies 

lke||L<^((0,T);L2(KiV)) = || || i2(]K]V) , ^^^^ 

e||V-ujLoo((o,T);L2(KiV)JV) < C (||-U°||2,2(i;]V) + e||VM°||L2(ig;]V)Jv) . 

Proof of Lemma 13. IL We multiply equation (jll) by 1*7 and we integrate by 
parts. Since all coefficients are real, taking the imaginary part yields 

d_ 
di 

Next we multiply (jlj by and we take the real part to get 

4/" (e^lVuA'^ + (c(-) +e^dJt,x)) \uf) dx = -2e^ f ^(t,x)\uJ'^ dx, 
dt J^N \ V V e / J J J^N at 

which yields (fTKll since is bounded in view of (jH). □ 

Our main result is the following homogenization theorem. 

Theorem 3.2 Assume ^ and ^ and that the initial data u'^ E H 



is 



n°(x)=V^„(^,rje2-^i;°(x), (16) 
with f° G H^{R^). Then the solution of ^\) can be written as 
u,{t, x) = e^^e2-^^„ v^{t, x) 



. e 



(17) 

+ e'-^e'^^^^m ( -, ^"^ ) vUt, x) + r,{t, x), 



■ e 

with ^ 

lim / / \r,{t,x)fdx = 0, (18) 

and (f„,fm) ^ C ([0, T]; L^(]R^))^ is the unique solution of the homogenized 
Schrddinger system 

.dv„ 



dt 



- div (AlVvn) + <^(t, x)v^ = in W X (0, T) 



z^-div(A;;,Vt;J+C(t,x)t;„ = x (0, T) (19) 

wjt = 0,x) = w°(x) inM^ 
L v^it = 0,x) = in M^, 
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with A* = ^Ve'VeXpiO^), for p = n,m, and 

dlJt, x) = dl^t, x) = \l dit, X, y)^^iy, e^)^^{y, 9^ dy. (20) 

Remark 3.3 Formula (H^, giving the coupling coefficient rf*^, is a ver- 
sion of the famous "Fermi golden rule " in quantum mechanics or solid state 
physics Jjjj, l^, \24]l - More precisely, the squared modulus of (i*^ is called 
the transition probability per unit time from state n to m and its formula is 
Fermi golden rule (see e.g. Chapter 6 in l^ ). The inverse tensor 
is called the effective mass of the particle corresponding to the wave function 
Vn J^, These effective coefficients not only depend on the chosen 

periodic crystal (characterized by the potential c{y)) but also on the energy 
level or Bloch band n,m of the particle, and on the quasi momentum 0""'"^. 
Effective mass theorems were already obtained in IS], chapter 4 of /F?]/ . 

However, the derivation of the coupled system 17^) and the justification 
of the Fermi golden rule is new to the best of our knowledge. 

Since 6^ and 6^ are not necessarily minimum points, the tensors A* 
and A*^ can be neither definite nor positive. Nevertheless, the homogenized 
problem ( TT^) is still well posed in C([0, T]; L^(M^))^ (by using semi-group 
theory J^), although its solution may not belong to L^((0, T); i7^(]R^))^. 

In general d*^ does not vanish, so that there is indeed a coupling between 
the two modes n and m. Nevertheless, if d{t,x,y) = d(t,x) in ^ and 
gn _ Qm ^ then, by orthogonality of the modes, we have rf*^ = 0. 

Remark 3.4 As already noticed in the scaling of (CP is not the usual 
semi- classical scaling for Schrodinger equation J^, \1 J]/ . \1 ^ . \14]l - The 

actual scaling of (Cp means that we are interested in much longer times than 
in the semi-classical limit. 

Remark 3.5 As already said in E.emark Wl\ all our results, including The- 
orem l/y.^ can be generalized if we replace the Laplacian by the more general 
operator diY{A{y)'V-) with a real symmetric, bounded, periodic and uniformly 
coercive matrix A{y). 

Remark 3.6 Theorem \S.S\ still holds true if the initial data is given by a 
combination of the two states 

= ^„ 9-) e'-'-^vl{x) + 9-) e'-'-^vl{x) 

instead of ( TT^) . Of course, it yields a non-zero initial data for Vm in ( TT^) . 
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Proof of Theorem 13. 2L This proof is in the spirit of our previous works 
[2], [SI- Define two sequences 



■->>n(9")t r,- 0"-x 

vl^{t,x) = u,{t,x)e ^ e ^'"^ - . 

Since = = by the a priori estimates of Lemma l3?T] we have, for 
p = n, m, 

Ike ||L°°((0,T);L2(]RiV)) + e|| Vt;f || l2((o,t) xR'V) < C, 

and applying the compactness of two-scale convergence (see Proposition 
I2.2p . up to a subsequence, for p = n,m, there exists a limit Wp(t,x,y) e 
((0,T) X M^;i7i(T^)) such that and eVwf two-scale converge to Wp 
and VyWp, respectively. Similarly, by definition of the initial data, v^{0,x) 
two-scale converges to ipn {y,9"')v^{x) and f™(0,x) two-scale converges to 
if ^ r and to tpn (y, v^{x) if 9"" = 9''. 
First step. We multiply ^ by the complex conjugate of 

e^(p{t,x,-)e'^^e^'''— 

where x, is a smooth test function defined on [0, T) x x T^, with 
compact support in (t, x) for fixed y. Integrating by parts this yields 

-le' [ u^Je-^'^'^dx-ie'' C [ v^^dtdx 
Jrn Jq J^n at 

+ [ / (eV + 2i7i9'')v';^ ■ (eV - 2z7r0")0' dt dx 
Jo Jr'^ 



+ / / (c' - A„(^") + e'd')v]^(f) dtdx = 0. 

Jo Jr'^ 

Passing to the two-scale limit yields the variational formulation of 



-(divj^ + 2^6*") (^(Vj; + 2m9'^)wr^ + c{y)wn = A„( 



)Wn in T'^. 



By the simplicity of A„(6'"), this implies that there exists a scalar function 
Vnit,x) e ((0,T) X M^) such that 

Wn{t, X, y) = x)tlJniy, 9"). (22) 
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Replacing n by m in the previous argument, a similar result holds true for 

Wm{t, X, y) = Vm{t, X)i>m{y, 6*™). 

Second step. We multiply JT]) by the complex conjugate of 



k=l 



k 



where (pit, x) is a smooth test function with compact support in [0, T) x 
M^, and Cniv) the solution of (fTTS|l . Integrating by parts, the resulting 
computation was already done in [B] in the absence of the oscillating potential 
d^{t,x). We briefly recall it: after some algebra, and using the summation 
convention for the repeated index k, we obtain 

r on nn 

Vu,-V^Jx= / (V + 2z7r— )(0t;,") ■ (V-2z7r— )^^ 



+e [ (V + 2^n^-){^v:) ■ (V - 2^n^-)C 
Jrn e oxk e 



e.|^< ■ (V - 2^n-)i;l 
oxk e 



+ [ (V + 2zvr^)(|^<)-e,V': (24) 
Jrn e oxk 



OXk 

+ / C'''(eV + 2^7rr)<-V- 



dxk 

A flrst simpliflcation arises from the deflnition of ipn which satisfles, for any 
smooth compactly supported test function $, 

p on nn i r 

/ (v + 2^7r-)V^^(V-2^7r-)$ + - / (c^ - A„(n)V^:$ = 0. (25) 
A second simpliflcation comes from the deflnition of 

/ (v + 2m-)e-(V-2^7r-)$ + - / (c^-A„(r))C'^^ = 

/ (V + 2^7r— )^^ ■ efc$ - e'^ / e^V-^ ■ (V - 2z7r— )$. 



(26) 
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Combining (|^ with the other terms of the variational formulation of 
we easily check that the first line of its right hand side cancels out because 
of (f25|) with $ = (^f", and the next three lines cancel out because of (f26|l 
with $ = ^^v^- We keep the three last terms of ([2^ which are bounded. 
Finally, JT]) multiplied by yields after simplification 



T 



v'J'W- — ■ ekipndtdx 

OXk 

<v4^ • (eV - 2t7ve'')Cn' dtdx 

OXk 

Jo Jrn OXk 
+ f f d'v'!'^, dt dx = 0. 







(27) 

We can pass to the two-scale limit in each term of (f27ll as was done in [H], 
except for the last one which is the only new and different term. The last 
line of (p7| is equal to 

/ a ( t, a;, — 1 - ( e ^ e ^ +e ^ e ^ 

JR^ ^ e/ 2 V 

u, + C(e)) dtdx 

= r [ d(t,x,-)Uv!''~'^ + vT) (i^^i-,9^mt,x) + Oie)) dtdx 
Jo Jr^ \ e/ Z \ e / 

where we introduced a new sequence ti^"-*^ defined, similarly to (f2T| . by 

„ .{2A,i{9")-Am(S'"))t „. (29"-9™')-x 

(t, x) = u,{t, x)e-' e'^^"- (28) 

Applying the same arguments as in the first step, this sequence t"^"-"^ is 
easily shown to two-scale converge to W2n-m(t,x,y) which satisfies 



(29) 



-(divy + 2i7r(2r - ^™)) (^(Vy + 2i7r(2r - e"'))w2n~ 

+C{y)w2n^m = (2A„(r) - XU0n)w2n-m m T^. 
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Because of the non-resonance assumption namely that 2Xn{0^) — 
is not equal to any eigenvalue Xp{29"' — 9"^), the spectral problem (f^ has no 
solution other than 0, which implies that W2n-m(t,x,y) = 0. The two-scale 
limit of (Ell) is thus 



I I \il^n?v^(t){t = Q)dxdy-i I [ [ {ipnl'^Vn^dt dx dy 

Rjv JjN Jq J^n Jjn at 

/ / i'nVn'^^ Ckipndtdxdy 

^M^ ^T^ ^^fc 



I I ipnVrN^-{Vy-2nxe^)cldtdxdy 
Jm.^ JT^ (J^k 



+ 



III Ct{^y + 2i7l0'')^nVn-V-P^dtdxdy 
Jo Jjn ^^k 

d(t, X, y)ipmVm4'n4^ dx dy = 0. 



1 '■^ 



2 Jo Jr'^ Jt'^ 

(30) 

To simplify (f3n|) we recall that J^^r \4'n\'^dy = 1, that d^^{t,x) is defined by 
(l2nll. and we introduce 



JjN V 

+^pnej ■ (Vy - 227rr)c' + ^nCk ■ (V, - 227re")cl (31) 
-C'(V, + 227rr)V'n ■ ej - Cl(V, + 2mr)^„ ■ efc)t/|/. 



Because of the Fredholm compatibility condition of equation (jlHl for the 
second derivative of ipm the matrix A*, defined by (f3T|) . is actually equal to 
8^V0V0A„(6'"). Finally JSHI) is equivalent to 



-i /" (pit = dx — i I I Vn-^dtdx— /" /" A* f „ ■ VV^fit rfx 
jR^f „io Jk^ dt Jq J^n 

d^^{t, x)vm(p dt dx = 0. 



T 




(32) 

A symmetric argument works for Vm (changing n in m in the test function 
\E'e). However, the initial condition is zero in the homogenized equation 
for Vm- Indeed, either 6*™ 7^ 9'^ and we already know that f™(0) two-scale 
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converges to 0, or 6'"* = 6'" and the orthogonality property 

'ipn'ipm dy = for m ^ n 



N 



implies that the coefficient in front of the test function 0(0) vanishes, which 
in the variational formulation implies that Wm(0) = 0. Therefore, (p^ and 
its counterpart for m instead of n yield a very weak form of the homogenized 
system (fT9| . Since A*^^A*^ are real symmetric matrices and d*^ = (fT9|l 
admits a unique solution in C ([0, T]; ^^(R^))'. By uniqueness of this solu- 
tion, the entire sequence two-scale converges weakly to t/'p (y, 6*^) fp(t, x) 
for p = n,m. 

It remains to prove the strong convergence (fTH|) . We compute 

lk<:(^)llL2(M'V) = lk€(^)lli2(]g]V) + \\i^n'^n(t)\\'j^2(^N-j + || "^m^m (^) II L2(IR'V) 
- 27^ / < it)lplvn it)dx-2n [ VT it)'^lvm (t) dx 

+ 27^/ 71 e^- 1 ^>„(t)^„t;„(t) rfx. 

(33) 

By the orthogonality property of the Bloch waves, the last integral in (|33|l 
converges to 0. By applying two-scale convergence, we can pass to the limit 
in the second line and in the last two terms of the first line of (jTH^ . For the 
remaining term we use Lemma f3. II which implies 

ll'"e('^)llL2(IRJV) = ||^e|lL2(KiV) ^ \'^nV^\\ll(^^N -^fN-^ = \\v°\\l2 



by the normalization condition of ipn- Thus we deduce 

lim llrf (t) II r2/iD.]V\ = l|t''^||r2/iD.]V\ — lltir ''"^'^ ""' ''■^'^"^ 



,,lL2(]g]V) — II llL2(]g]V) - ||t^nl'';|lL2(]g]V) - ||tVmVt;||^2 

which is precisely because of the conservation of energy of the homogenized 
system (fT9|) . i.e. 

ll'^n(^)llL2(RJV) + ll'"m(^)llL2(R'V) = II'^°IIl2 



Since lim^^o lke(i)||^2(]j]V) = 0, the Lebesgue dominated convergence theorem 
yields □ 
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Remark 3.7 Recall that the function (kiv) is the solution of / fI3|) . unique 
up to the addition of a multiple of ipn- This multiple may depend on (t,x) 
and therefore the test function as well as the homogenized system could 
depend on the choice of this additive term. Actually the homogenized system 
depends on (k only through the homogenized tensor A'^, defined by ^Sl\) . If we 
replace Ck{y) by Ck{y) + Ck{t, x)il)n{y), an easy calculation shows that all terms 
Ck cancel out because of the Fredholm alternative for (k, i-e. the right-hand 
side of is orthogonal to ipn- Thus, the homogenized system is uniquely 
defined whatever the choice of the additive constant in Ck{y)- 

Remark 3.8 A formal two-scale asymptotic expansion (in the spirit of ^) 
of the solution of (CJj would give 



As usual in periodic homogenization, this expansion suggests the choice of the 
test function in the proof of Theorem \S.^ Another possible interpretation 
of is as follows. The large terms in the variational formulation of 
(CP cancel out because of the equation satisfied by ipn. However, new terms 
of order e^^ appear because of the first order derivatives of ipn- They are 
compensated in turn by the second order derivatives of the corrector C^. 

Remark 3.9 Part (i) of assumption ^ states that the eigenvalues \n{0^) 
and \m{0"^) are simple. This hypothesis is crucial in order to be able to 
differentiate the spectral cell problem with respect to 9. If one of these eigen- 
values is not simple then, as is well known, it is not anymore differentiable, 
but merely directionally differentiable (which is not enough for our purpose). 
So, we do not know how to generalize Theorem \S.S\ in the case of multiple 
eigenvalues. There is one notable exception when one eigenvalue is of mul- 
tiplicity, say p > 1, and there exists locally a labelling of the eigenvalues and 
eigenvectors in p smooth branches. Note that it is a very strong assumption, 
which is rarely meet in practice. Then, using an argument of one can 
generalize Theorem \S.2 and obtain a limit system similar to with as 

many equations as the repeated multiplicities of the eigenvalues A„(6'") and 
Am(^™'); and coupled only by zero-order terms. 



Ue(t, x) ~ e 




) 
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Remark 3.10 Part (ii) of assumption ^ states that the group velocities 
vanish, VeXn^O"') = VeAm(6'™) =0. If it is not the case, then it induces 
a large drift of order and the homogenized system ^^) can be obtained 
only in a moving frame of reference, following this large drift (see for 
more details). Therefore, if 'VeXniO'^) 7^ VeiAm(6'™), one can not generalize 
Theorem LV.^ since both initial and target states move with large different 
speeds, so no coupling is possible in the limit as e goes to zero. In the case 
V6)A„(^"') = VeXmiO"^) ^ it is technically possible to generalize Theorem 
\S.S\ following the argument of J^, but this result would not make much sense 
since it would assume that the exterior potential de{t,x) move with the same 
velocity, or at least is macroscopically constant, which is usually not the case 
in physical applications. 



4 Electromagnetic potential 

Instead of (QJ) we now consider a Schrodinger equation with an exterior elec- 
tromagnetic field 

i-^ - (div + ieAe) (V + ieA^) u, + e'^c j u, = in x (0, T) 
u,{t = 0,x) =u%x) inM^, 

(34) 

where Ae{t,x) is the electromagnetic vector potential, i.e. a function from 
M+ X into 0, PI, [21. The electric field E and magnetic field B are 
recovered by 

dAe 

E{t,x) = — -^{t,x) and B{t,x) = cuT\Ae{t,x). 

For an electromagnetic wave, the vector potential is assumed to be given by 

A,{t, x)=n [e' ^ e^'^ : j a (t, x, - j , (35) 

where TZ denotes the real part and a{t,x,y) is a bounded smooth function 
from M"*" X x into M^. As before, c{y) is a bounded function from 
into M, the initial data m° belongs to if^(]R^), and the conclusion of Lemma 
I3.1l still holds true: there exists a unique solution of (fHij) in C ([0, T]; iJ^(]R^)) 
which is uniformly bounded in ((0,T) x R^), independently of e. Theo- 
rem [HSl can be generalized as follows. 
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Theorem 4.1 Assume ^ and ^ and that the initial data vP^ G H^{M.^) is 



with e H\R^). The solution of (M) can be written as 



u,{t, x) = e^^e^^'^^V^^ (-, ] Vnit, x) 



e e 



''■ipm [-,0"')Vm{t,x) +r,{t,x), 



with 



lim / / |re()f:,x)| dx = 0, 

,2 



and {vn,Vm) G C ([0, T]; L^(]R^)) is the unique solution of the homogenized 
Schrddinger system 



- div {A*^Vvn) + dl^{t, x)v^ = Q in x (0, T) 

dv 

z-^ - div {A;,Vu^) + dl^^t, x)vn = ^ in x (0, T) 



dt 

Vnit = 0, a;) = f "(a;) 
^ v^{t = 0,x) = 

with A* = e'^piP'^) , for p = n,m, and 



(36) 



m 



m 



JN 



/ i^Ay. on^it, X, y) ■ (V + 2inen iJm{y, on dy. 

(37) 

Remark 4.2 /n general d*^^ does not vanish, even if a{t,x,y) is a constant 
vector and ^" = 6"^, so that there is indeed a coupling between the two modes 
n and m. 



Proof of Theorem 14. IL The proof is very similar to that of Theorem 
13 .21 The first step is identical, and in the second step we choose the same 
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test function defined by (f^ . The higher order term in the variational 
formulation is 

/ (V + ieA,) ■ (V - ieA,) '^Jx = / Vu, ■ V^, 

+ie / ■ V"*, - If.A ■ Vm,) rfx 



/ lAefu^'^^dx. 

(38) 

The first term in the right hand side of (j38| is exactly the previous term (f24| . 
The last one goes to zero, while the second one is the only new term which 
yields a non-zero limit. Indeed, integrating by parts in this term gives 



dx 



ie / {u,Ae ■ V^, - ■ Vm,) dx = ie / u, {2A^ ■ V^, + ^,divA) 

= i/2 I -"^ + v'^^ (2a' ■ {Vy - 2i7r^")^; + ^^div^a^) dx 
+i/2 I 2m (^;2"-'»(^'» _ . a' + - 9"^) ■ a') H^^dx + 0{e). 

(39) 

Recalling that the two-scale limit of ff" is 0, the limit of (j39|) is 

i/2 I I vjp^rn (2a • {Vy - 2i7r^")V^„ + ^„divj,a + 2i7rV^„(^" - ^™) ■ a) rfx dy 

which yields formula P7|l for the coupling coefficient d*^. The rest of the 
proof is identical to that of Theorem 13.21 □ 

5 The resonant case 

In this section we come back to the original Schrodinger equation JT]) but we 
change assumption ^ by assuming that there is a single resonance between 
the initial data and the target state, namely 

{Hi) there exists / > 1 such that Xi{29'^ - ^™) = 2A„(0") - A„(0'"), 
(iv) foranyp>l, Ap(3r - 20^") ^ 3A„(r) - 2A„(0'"). 

(40) 
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We keep assumption ^ that we extend to the new eigenvalue A/ for the 
Bloch parameter 6^ = 26^ — 9^, i.e. 



for p = n,m,l 



(i) \{d^) is a simple eigenvalue, 
{ii) 6*^ is a critical point of \{6) i.e., Ve\{0^) = 0. 

(41) 

With these new assumptions we generalize Theorem 13.21 hv obtaining a limit 
system coupling three possible states instead of just two. 

Theorem 5.1 Assume and and that the initial data G H^(M.'^) 
is 

with e H^{R^). The solution of ^\) can be written as 



ipi (^,^') vi{t,x) +r,{t,x), 



with 

rT 



lim / / \r^{t,x)f dx = 0, 
Jq J^n 

and {vn,Vm,vi) G C ([0, T]; L^(M^))'^ is the unique solution of the homoge- 
nized Schrddinger system 

dv 

t-^ - div {AlVvn) + dl^{t, x) Vm + d^t, x)vi = in x (0, T) 
dv 

i-^ - div {Aiyvm) + d*^St. x)vn = Q in x (0, T) 

i-^ - div {A*iVvi) + d*i^{t, x) w„ = in x (0, T) 

= 0,x) = mM^ 

f™(t = 0,x) = mM^ 

vi{t = 0,x) = in M^, 
with A*p = eV eXpiO''') , for p = n,m, I, and 

1 
2 



x) = 4„(i, a;) = - / d{t, X, y)ij^iy, e^)^p{y, 9^) dy (42) 
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for p = m,l. 

Remark 5.2 More generally, there could be multiple resonances between the 
initial and target state. Let ko > 1 be the order of the resonance. Under 
a suitable generalization of assumption UOj ), all modes of momentum {k + 
1)^ - kO"" and energy {k + 1)A„(^") - A;A„(^™) are coupled for -1 < k < 
ko. Theorem, \ 5.1\ can be generalized to obtain an homogenized system for 
{vm,Vn,V2n-m, ■■■,V(^ko+i)n-kom) in wMch thc coupUng matrix d* is hermitian 
of size ko + 2 with the following sparse structure 
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Proof of Theorem 15. IL The only modification with respect to the proof of 
Theorem 1221 is the fact that the sequence ti^"-"*^ defined by now admits 
a non-zero two-scale limit ipi{ii^6^)vi{t,x) because the spectral cell problem 
(f29|) has a non-trivial solution t/)/, as a consequence of part {iii) of assumption 
(|in|) . No other states appear because of part {iv) in (|^ . The rest of the 
proof is similar to that of Theorem 13.21 and we safely leave it to the reader. 
□ 
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